Dipolar bogolons: from superfiuids to Pfaffians 
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We study the structure of Bogoliubov quasiparticles, 'bogolons,' the fermionic excitations of paired 
superfiuids that arise from fermion (BCS) pairing, including neutral superfiuids, superconductors, 
and paired quantum Hall states. The naive construction of a stationary quasiparticle in which 
the deformation of the pair field is neglected leads to a contradiction: it carries a net electrical 
current even though it does not move. However, treating the pair field self-consistently resolves this 
■ problem: In a neutral superfluid, a dipolar current pattern is associated with the quasiparticle for 

D " which the total current vanishes. When Maxwell electrodynamics is included, as appropriate to a 

Ph \ superconductor, this pattern is confined over a penetration depth. For paired quantum Hall states 

of composite fermions, the Maxwell term is replaced by a Chern-Simons term, which leads to a 
dipolar charge distribution and consequently to a dipolar current pattern. 
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Introduction: Paired superfiuids are among the most 
ubiquitous of the many ordered phases of interacting 
fermions in two and three dimensions. In condensed mat- 
ter settings they include both neutral superfiuids such as 
3 He, charged superconductors and now also paired quan- 
tum Hall liquids such as the Moore-Read or Pfaffian state 
that is believe to underlie the quantized Hall plateau at 
filling factor v = 5/2. In all cases paired superfiuids 
exhibit two distinct excitations that dominate much of 
their physics: vortices and Bogoliubov quasiparticles or 
'bogolons'. The former are a generic consequence of su- 
perfluidity but the latter are a particular signature of 
pairing — they involve breaking apart a Cooper pair into 
its fermionic constituents. 

The structure of vortices is well understood: They are 
the topological solitons of a complex scalar order param- 
eter in the Landau-Ginzburg description of a superfluid. 
In a superconductor, additional coupling to a Maxwell 
gauge field results in an associated quantum of flux, while 
for a quantum Hall liquid, coupling to a Chern-Simons 
gauge field associates a quantized charge with each vor- 
tex. The structure of bogolons is less well understood as 
they are, by comparison, much more quantum mechan- 
ical particles. We will address that gap by providing a 
theoretical analysis of their structure for all three exam- 
ples alluded to above. For superfiuids and superconduc- 
tors we will be able to recover the heuristic description 
advanced by Kivelson and Rokhsar [l| . For paired quan- 
tum Hall liquids our results are new and add to a recent 
burst of interest in the properties of bogolons @, Q , in- 
cluding work by four of the present authors 0] ■ 

In the weak pairing (BCS) limit the momentum (or 
Bloch) eigenstates of the bogolon exhibit the well known 
dispersion relation sketched in Figure la, with a charac- 
teristic minimum at the underlying Fermi surface. In 
terms of these, one can make a localized wave-packet 



state with a spatial extent large compared to the coher- 
ence length, £, and a well defined momentum. Unlike a 
wavepacket in the normal state, this bogolon wavepacket 
has a group velocity which is different than the Fermi 
velocity Vp, and which vanishes on the Fermi surface. 
It has spin 1/2, but its (average) charge is smaller than 
the charge of electron e. Both quantities vanish as the 
momentum of quasiparticles p approach the Fermi mo- 
mentum pp. On the other hand, since the wavepacket 
has a net momentum, it carries a net current [H[ equal 
to evp . This indicates that our construction of a local- 
ized bogolon is fundamentally inadequate. The problem 
becomes especially clear in the limit p = pp, where the 
group velocity of the wavepacket is zero. In this case the 
current density is finite inside the wavepacket and zero 
outside of it. The resolution of this puzzle will lead us to 
a bogolon structure that involves an algebraically falling, 
dipolar, return current flow via the condensate for neu- 
tral superfiuids, a version of this screened on the scale 
of the London length for superconductors, and a version 
exhibiting a charge dipole as well as a locally dipolar 
backflow for two dimensional quantum Hall fluids. Alto- 
gether, bogolons are fairly complicated objects! 

Bogolon wavepacket: We begin with the mean-field 
BCS Hamiltonian for a neutral fully gapped paired su- 
perfluid which also serves to fix notation, 



h bcs = Yl &<L<** + E [ A k4t c -H + h - c -] w 
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where £ k = — fi, A k is the gap function and we work in 
units where H = e/c = 1. It is a simple matter to diago- 
nalize, H BCS = £ k , s £k7k, s 7k,si with ^k = + |A k | 2 
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FIG. 1: (a) Quasiparticlc dispersion, (b) Current flow 
around neutral supcrfluid bogolon (window size ~ 2A.) 



by means of a Bogoliubov transformation, 

■A 



7k,t = "kCk.t - u kcl k ^,7k4 = Uk4;,t + U kC-k,i, 



in 
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The (T = 0) BCS ground state is then the state annihi- 
lated by all the 7 ks , \Q) = TIk>o ( u k + WkC k t c -k|) I )- 
Singlc-quasiparticle states with momentum k are given 
by |ks) = 7 k J^) and it is readily verified that their en- 
ergy _Bk is minimal at |k| = kp. We will work in d = 2 as 
that naturally includes the case of the paired QH state, 
but the results are readily generalized to d = 3. 

A quasiparticlc wavepacket with average momentum 
fiko = hkpko, spin s and spatial extent ~ A is obtained 
by superposing the states |ks) with momenta near ko 
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In order that the energy uncertainty of the wavepacket be 
smaller than its average energy, we need to choose A 3> 
£ = as can be deduced from the low lying dispersion 



[v F ko-(k-ka)]- 

2An 



where An = |Ap 



relation E(k) « A + 
and = kp /m is the Fermi velocity. 

Our primary concern is the structure of quasiparticlc 
wavepackets centered at momenta close to pp, so that 
their group velocity is much smaller than vp. Clearly 
the packet has vanishing group velocity at p = pp. How- 
ever, a tedious but straightforward computation of the 
expectation value of the quasiparticle current operator 



iqp 

Jq 



Z^k,s in k+3 s 



Ck-s s in the state yields 



5 qp \- T . = v F k e 



(4) 



We are thus presented with a contradiction: a stationary 
quasiparticle wavepacket is associated with a current that 
has nonzero divergence — violating the continuity equa- 
tion. 



A first step in resolving this puzzle is to observe that 
we have taken a slippery step in passing from momentum 
space to real space. In real space, the wavepacket state 
([3]) is now inhomogeneous and hence a homogenous "pair 
potential" A no longer yields a self-consistent mean field 
theory of the wavepacket 0- It is possible to prove that 
any state that satisfies the self-consistency conditions re- 
spects the equation of continuity. Recomputing the pair 
potential in the wavepacket state and then iterating the 
construction of the wavepacket and the computation of 
the pair potential should yield a state that does obey 
current conservation Q. In the supplementary material, 
we show that the first iteration of this process produces a 
change in the pair potential that already partially cancels 
the quasiparticle current. 

However implementing this approach requires detailed 
numerical work. Instead, we construct an effective action 
which correctly treats the low-energy, long-wave-length 
physics in the weak coupling limit, An, <C Ep = k F /2m. 
While portions of this work may be reconstructed from 
existing literature, in particular the 'conserving approx- 
imations' to superconducting response, to our 
knowledge an explicit quantitative treatment of a bo- 
golon wavepacket has not been previously presented. 
Neutral superfluids: As we are interested in a 
wavepacket constructed from momenta very close to the 
(parent) Fermi surface, it is sufficient that we work with 
the effective dynamics for this set of degrees of freedom. 
Formally, we begin with a Hubbard-Stratonovich (HS) 
decoupling of an attractive four-fermion interaction in 
the particle-particle channel, and integrate out fermions 
above a cutoff thus generating an effective action for the 
HS field, A(r,f) = A e ie ( r '*). As we are in the broken- 
symmetry phase, fluctuations of the amplitude can be 
neglected. The result is an effective theory of dynamical 
fermions coupled to a dynamical phase field 0(r,t) [H |. 

To be explicit, we consider the case of s-wave pairing, 
where the most important terms in this (well known) 
theory are represented by the action S = J dtd 2 r {C^ + 
£ p + Cg), with 
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4(r,t) , (5) 



C p = -Ane" 



tyj(r,t)^(r,f)+h.. 
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where xo is the static compressibility (equal to the den- 
sity of states at the Fermi surface), and n s is the super- 
fluid density. At T = 0, n s = p, the total electronic den- 
sity, for the Galilean invariant systems considered here. 

Note that the conserved charge is no longer carried 
solely by the fermions, but also by the supcrfluid compo- 
nent via twists in the order parameter. A straightforward 
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application of the Noether procedure allows us to write, 
for the density and current 



P 



Zm 



where pf* = £ s and j qp = £ s ^M*I>1 V^ s 
S we then obtain the equations of motion 



d,p" = -Vf 



Br, 



(8) 
From 

(9) 
(10) 



where B p = 2zAq 



e t6 ip^ip^ — e ie V , 4-V'tJ i s ^he term that 
couples the quasiparticles and the superfiuid. From ©, 
© and (flO]), it is evident that d t p + V • j = 0, i.e. the 
properly defined density and current obey the continuity 
equation; it is equally clear that the quasiparticle density 
is not independently conserved. 

Let us now specialize to the treatment of a stationary 
bogolon wavepacket in the approximation where we ig- 
nore the quantum fluctuations of 9. This implies that the 
LHS of Eqns © and (jTOj) vanish so that 
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Thus, in the wave-packet state for which (j^ 15 )* is given 
by Eq. [H the resulting phase texture is 



i{q- ko) 
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which permits us to write for the total current 



(jq) 



VF 



q 2 k - (q- k )q 



(12) 
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Eq. (jTB")) corresponds to a real space current (j(r))$ 

z x V^aW, where (p x (r) = 2nv F {%0 ^ 2 (l - e^l* 

The flow pattern is solenoidal (clearly V • (j)* = 0), 
and decays as r~ 2 far from the center of the wavepacket. 
Corrections to this expression at short distances are non- 
universal, and are beyond the reach of the field-theory 
approach. Finally, we note that at finite quasiparticle 
concentration p qp , the long-range nature of the distri- 
bution of current density in a quasiparticle wave packet 
leads to the conventional expression j qp = evFP qp for 
the quasiparticle contribution to the current density, in 
agreement with the Boltzmann approach [l5j applicable 
in this limit. 

Superconductors: We now turn to the case of a 
charged superfiuid which is minimally coupled to a fluc- 
tuating U(l) Maxwell gauge field - i.e., the super- 
conductor with dynamical electromagnetism. The effec- 
tive action is obtained from that of the neutral super- 
fluid by converting the derivatives to covariant deriva- 
dn — iAu, where the dynamics 



tives: -> — ^ - 

of A^ are described by /^Maxwell 



jF^F^ in which 



Fpv = d^Ay — d v A^ is the Maxwell field strength. From 
S + ^Maxwell, we find the equations of motion for the 
quasiparticle and superfiuid currents 

P = P qP - Xo(d t - 2A ), j = j clp + ^(Ve - 2A) 



d tP q 



B„ 



Xo d t (d t 9-2A ) = -^V- (ye - 2A) + B p (14) 
2m 

supplemented by Maxwell's equations 

V • E = 4tt(p — p), V-B = 0, 
V x B = 47rj + d t E, VxE = -d t B. 



(15) 



In and (fT5|) the quasiparticle current and density 
take their gauge-invariant forms, p qp = ^ s iplifi s and 
j qp = E s ^Im[^D^], and E = -d t A - VA a and 
B = V x A are the electric and magnetic fields in the 
quasiparticle state; in writing the Poisson equation we 
have assumed the existence of a neutralizing positive 
background p — (J2 S ^\i>s)n hi the BCS ground state. 

The first comment to be made here is that now even 
extended bogolon states of definite momentum do not 
carry current. This basically reflects the Meissner ef- 
fect. Specifically, the uniform quasiparticle contribu- 
tion to the current is exactly cancelled by a superfiuid 
backflow, which in unitary gauge 6 = 0, corresponds to 
^A = (j qp ) oc Vpko [Iff. The correct bogolon state 
carries no current; they are neutral particles. 

Still in unitary gauge, let us turn to the construction 
of the wavepacket. For static wavepackets we find that 
the third equation of (H) yields • A = (B p ) = (V • 
j qp ), so that as before the total current j = j qp — —A 
is conserved. Using this, we rewrite the third Maxwell 
equation as 

[-V 2 + A^ 2 ] A = At: (j qp - A 2 V(V • j qp )) (16) 

where we have defined the penetration depth via A^ 2 = 
™" and the expectation value is taken in the naive wave 
packet state with A^ = A = 0. Using either the expec- 
tation value of (Bp) computed in the superfiuid case or 
the form of j qp , we may solve (|16[) by Fourier analysis: 



(jq)* = v F 



q 2 k - (q- k Q )q 

i 2 + K 2 



(17) 



which coincides with ([131 in the limit in which the cou- 
pling to electromagnetism vanishes (when Al — > oo). It 
is easy to see that the power-law asymptotics of the su- 
perfiuid case are replaced by exponential behavior at long 
distances, j(r) ~ e ~ r / x L for r > A and Al. This reflects 
the fact that superconductors screen magnetic fields and 
thus the current pattern is confined to within a penetra- 
tion depth of the center of the bogolon (Note that 
the short-distance behavior of the wavepacket is qualita- 
tively different depending on whether the superconductor 
is Type I (or weakly Type II) in which case A l completely 
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characterizes the current distribution, or strongly Type 
II, with £o "C A <C Xl, in which case the bogolon resem- 
bles that in a neutral superfluid for r -C A^.) 
Paired QH States: Our final example is the case of 
a bogolon in a paired QH state of composite fermions 
(CFs). Here, we start with fermions moving in a static 
uniform background field A (where V x A = B), and per- 
form a 'flux attachment' by means of a statistical gauge 
field a whose dynamics are governed by a Chern-Simons 
(CS) term, Ccs = j^e^ up apd u a p , with 3> the quantum 
of flux. Qualitatively, the role of the CS gauge field is 
to attach two quanta of magnetic flux to each electron 
to convert it into a CF, which sees zero net flux at half- 
filling, i.e. we have B = 2$op~. In this case, we replace 
dp — > Dp = dp — i(a + A)p, and change the currents 
and densities accordingly. Although more properly we 
should consider the example of spinless fermions and a 
p-wave superconductor, the distinction is unimportant as 
we are primarily interested in the interplay of the CS elec- 
trodynamics and charge conservation, neither of which 
depends essentially on the pairing symmetry. The equa- 
tions of motion now follow as a result of S + Scs : the 
'matter' equations arc similar to the previous example, 

P = P qp - XoldtO - 2(a + A )}, (18) 

71 

J = J qp + ^ [V0 - 2(a + A)] , dtp** = - V • j«» + B p 

71 

Xo d t [d t - 2(a + A )\ = -i-V • [V0 - 2(a + A)] + Bp, 

2m 

but the Maxwell equations are replaced by the CS equa- 
tions, which are pure constraints: 



b ee V x a = -2$ { P w - X o[d t 9 - 2(a 

n 



e ee — d t a — Va = 2$ z x 



{j qp 



2m 



Ao)]} (19) 
[V0-2(a + A)]}. 



Note that now, A is not a dynamical field, but rather 
represents the background magnetic field, - 



2*o 



P, 



where p is the mean density, (pf)o- in the groundstatc, 
and, for the present, we will set the external potential 
Aq = [lij . Proceeding to the wave packet state, and 
specializing to unitary gauge and to static configura- 
tions as in the previous examples we use the identity 
V 2 V = z x V[V x V] + V[V • V] (valid in d = 2) to write 

[-V 2 + A" 2 ] (a + A) = 8 Xo ^ <j qp - A 2 CS V(V • j<*)) 
where A™ 2 = 8x °^ " 3 . The solution for the total field 

^s m 

a + A is similar to (fT7|> ; observe that the flux is expo- 
nentially screened over a distance A C s- Using \o ~ ^ 
as appropriate to a 2D Fermi surface, and the QH re- 
lation p = l/4ir£ B for filling factor v — 1/2 where £b 
is the magnetic length, we find A cs ~ \Ib (p/ns) 1 ^ 2 , so 
for a Galilean invariant system at T = where n s = p, 
we find Acs ~ \^b- Thus, the characteristic size of a 
bogolon wave packet is of order the magnetic length. 
A striking difference from the normal superconductor 
is that the second CS equation forces the existence of 



an electric field, which leads to a deviation of charge 
density from the background. The simplest estimate is 
Sp ~ -kr^^e~ 2r ^ B \ while this is not the exact form, 

t B r 

the important point is that there is necessarily a dipo- 
lar charge distribution oriented perpendicular to fc , with 
separation ~ fg, accompanying the screened dipolar cur- 
rent pattern. Upon inclusion of the long-range Coulomb 
interaction (ignored so far) both current and charge 
densities acquire power-law tails similar to those in [20( . 

In the QH case, the bogolon has a natural interpreta- 
tion as the descendant of the CF in the paired phase. Sev- 
eral authors, including one of us 2ll424|. have observed 
that upon projection to the lowest Landau level the CF in 
the compressible phase goes from being a charged particle 
to a neutral particle with a dipole moment proportional 
to its speed and perpendicular to its direction of propaga- 
tion. The argument for a dipolar charge distribution for 
the bogolon presented here - the application of CS elec- 
trodynamics to a paired superfluid - is rather different 
from projection to a reduced Hilbert space, and the con- 
nection between the two cases is an intriguing question 
that we hope to address in the future. 
Concluding Remarks: In this paper, we have given a 
consistent microscopic description of bogolon wave pack- 
ets in three broad classes of paired fermion states: super- 
fluids, superconductors, and paired composite Fermi liq- 
uids with CS electrodynamics. In all cases, the quasipar- 
ticle is associated with a decidedly nontrivial current flow 
pattern carried in part by the condensate, and manifestly 
obeys global and/or local conservation laws as appropri- 
ate. Although for pedagogical simplicity we focused on 
the case of stationary wave packets, this restriction is 
merely a matter of convenience: suitably boosted current 
configurations are associated with bogolons in motion. 

Our results are valid in the limit p^P^ 2 <C 1 when the 
concentration of quasiparticlc wavepackets is small , or in 
other words when the distance between quasiparticles is 
much larger than their size. In the opposite limit of high 
quasiparticle concentration where quasiparticles overlap, 
the system can be studied using the kinetic equation ap- 
proach 15(. In this formalism, the Boltzmann equation 



for the distribution function of quasiparticles is sup- 
plemented by equations of motion for the electrodynamic 
fields and continuity equations expressing charge conser- 
vation. The current and charge densities take the form 



j = epv s + e I d I k — rik 



p = e 



J d d k[uln k + vl(l-n- k )}, 



(20) 



where v s = ^ (V0 — 2A) is the superfluid velocity. The 
situation is similar to the microscopic scenario discussed 
here: in order to describe the distribution of v s , an addi- 
tional variable included in the kinetic theory (compared 
to the case of the normal metal), charge conservation 
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must be treated as an independent equation, rather than 
following directly from the equations of motion. 
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SUPPLEMENTARY MATERIAL 



Bogoliubov-de Gennes Equations and the BCS Ground State 



We begin with the reduced BCS Hamiltonian, 



#bcs = Y £ kC L c ks + Y AkC kt c -i4 + A k c -k4.c kt 



(21) 



k 

«=t,4. 



where £k = ^7 — £t, and Ao is the pair potential. As is standard, we define Bogoliubov operators via 

7k.t = "kCk,t - vuc f _ Kl ,j^i = «kc k ,t + u kC-k.i 



(22) 



Note that for the quasiparticle operators we have assigned spin by noting that removing a spin up and adding a spin 
down both lead to a state with a net spin of —1/2. As usual, we determine «k, ^k by requiring the Hamiltonian to be 
diagonal in the quasiparticle operators, which gives 



1 



1 



(23) 
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where £ k = + |A k | 2 , and 



k 



Since |iik| + l^k| — lj we may invert the Bogoliubov transformations as 

Ck,t = Wk7k,t + «k7 k ,i> c -k,l = -«k7k,t + u k7k4 



(24) 



(25) 



The eigenvalues of the BdG hamiltonian come in pairs with positive and negative energies. The T = BCS ground 
state is obtained by filling all the negative-energy quasiparticle levels 



1^ = 11(^+^4^) |o) 



k>0 



(26) 



and is annihilated by all the Bogoliubov operators, 7ks |r2) = for all k. 

The Wavepacket State 

The quasiparticle dispersion has a minimum at k = hp- An excited state with one quasiparticle with momentum 
k and energy and spin s is given by 

|k S )= 7k JO) (27) 
A low-energy wavepacket solution with average momentum kg and spatial extent ~ A can also be constructed, as 



I k ,s/ 



d rf fce -iA 2 (k-k ) 2 | ks 



d d fce -iA 2 (k-k ) 2 7 tj O) 



(28) 



with |ko| ~ hp, so that the energy is small. For ||ko| — "C Ao, the energy is approximately E(ko) w Ao 
c ^ v F {\]ia\-kF) wnere Cq j g a num k er f order 1. 



Current and Density in the Wavepacket State 



The current operator is given by j q = ^ ks £c k+as c k _a s , and 



«.s\k\<,s) = (^)7 ^'j ^fc" e ^ A2 [ (k '- k0)2 + (k "- k0)2 ](k' S |j q |k".s) 



(29) 



We may now reduce the current matrix element in the excited state to a ground state average, using the fact that the 
Bogoliubov operators annihilate the BCS ground state |f2): 

(k' S |j q |k" S ) = <fi|7k' s j q 7k|tt) 



(fi| 7k , s ( 



Jq:7k 



Jq:7k 



+ 7k»Jqj \ n ) 

}|fi> + <fi|j q |fi>*(k'-k") 



(30) 



The first term may be simplified as 



{7k' t , [j q , 7k'' t ] } = J2 — {"k'Ck't ~ «k'C f _ kU 

k,s 



c k+ a s c k _a s , ii k „c k „ it -U k ,c_k''4 



} 



E — { u k'C k , t - Wk'cL k , 4 , w k »c k+ a t <5 (k" - k + |) + u k „Ck_a 4 5 (-k" - k - |) } 

k 

* — ' m I 



k 

k' + k" 

2m 



2/ V 2 
(M k 'M k » + «k'W k ») S(k' - k" - q) 



uvuL,6 k'-k-^Uk"-k+^- wvt„5 ( -k' - k + ^ ) 5 I -k" - k 



(31) 
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and, similarly 



f '11 k' + k" 
[ik'i, [jq, 7k»|J } = ( M k'"k'' + «k'«k") <5(k' - k" + q) 



(32) 



while the second term can be made to vanish by choosing a real wavefunction for the ground state (check?) 
For concreteness, let us add a spin up quasiparticle. Then, the current in the wavepacket is 



« t |j q K,T> ={^) J d ' 



X 



v d fce -A 2 (k-k ) 2 + ^l 



k 



(33) 



For the s-wave case, we may take (in the first iteration, i.e. at the mean-field level before enforcing self-consistency) 
the ms and vs real, in which case the purely q = contributions to B p and j both vanish. For specificity, we work in 
two dimensions and choose a stationary wavepacket with k = ipi, which corresponds to a quasiparticle at rest[25|. 
If wc define k = k + k, we have 



l«k±a 



j_ J v F (X ± f) 
V2 1 1 + 2A 



0(fc 2 ),kk±a| 



1 f v F (h±*r) 



2A 



0(fc 2 



which, when substituted into the expression of the current (with k shifted as above) yields 

2 



Uq)* = C 



X 



2T„-\ 2 k 2 1 



d z k 



fcpX + k 



j ( v F (h + \)\\ »F (K -t)\ if «F (k + t ) } f «F (*. - f) 



For the density in the wave packet state, a similar calculation yields 



(*k ,t|pq|*k ,t>=^q,0 



A 



(34) 



(35) 



(36) 



Backflow Source Term 

Th source term to compute the backflow current is given by the expectation value of the operator 
B p (v) = 2tA (v}(r)Vl(r) - 1>&)M*)} 



2^oE {4,A',^ (k - k,) - r - - .'v- k k ' ■•} 

kk' 

2*A ^ { (« k7 t it + < 7k4 ) (-<, 7k ',t + «k'7 k%i ) ^ (k - k,) ' r 



-«k'7^,t + w k'7k',i) (<7k. t + AJ e^ l(k " k,) - r 



e-^ k - k ' ) r } 



= 2^A ^ { («k7k,t + w k7k,|) (-<-7k',t + "k'7 k .i) - ("At + <7k,+) («kOV,t + } e l(k k ' } - r 

kk' 

= 2iA ^ |(«k'Uk - v^,u k ) 7 kjt 7v,t + « w k' - w k Wk') 7k,i7k',i 

kk' 

+ K« k ' + tfctv) 7k, t 7k'4 + (-«' - Wk«t) 7k,47k',f} e 4 ( k - k ')- r (37) 
where in going from the third to the fourth line we have interchanged k, k' in the second term. 
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Proceeding as before, we compute expectation values of the source term in the wavepacket state, and find that only 
the first term has a generic contribution; the second term contributes only when k' = k", while the last two vanish: 

d 



(B p (r))« = (-^=j J d d k'd d k , 'e-^ 2 ^'~^ 2+ ^''-^(k' s \B p (r)\k' , s ) 
= 2tA (-j=y I d d k'd d k"e-^ k '-^ 2+(k "- k ^ [( 



u k „v k > ~ u k //Uk'J 
+<J(k' - k") (v* k ,u w „ - ut,v k „)] e^ k '- k "^ (38) 



We may rewrite this in momentum space, (B p (r))^ = J j d h (B Pq )^e lq r , whence 



(# Pq )* 



2iA n 



d d fce -^(k-k„) 



u k +^k-f 



(39) 



Using the expressions Q34[) for it and v and for simplicity [26j working in d = 2 we have, explicitly 



(B Pq )* = 2iA e~ 



* 2 <, 2 / A 




vf ( fcx + % 



2An 




t>f (X - f ) 1 if ( 
1 2A [ 2 I 1 2A 




vf (k - f ) 
1+ 2A 



iq x v F e 4 



(40) 



Gap Function and Phase Texturing in the Wavepacket State 



We now show that the expectation value of the gap operator acquires a phase winding in position space when it is 
computed in the wavepacket state built from the uniform BdG solutions above. We have 

A*(r) = V>j(r)V{(r) 

— C k,t L -k'4 e 



kk' 



kk' 

{- w k' w k7k, t 7k',t + u k u k'7k,;7k'4 + w kUk'7 k , t 7 k '4 _ «'7k,.L7k',t} e l(k - k,) - r 



(41) 



kk' 



In the wavepacket state, we proceed as in the previous section and find as before that the latter two terms vanish, the 
second is a purely q = contribution, while the first term is 'generic'. Going to momentum space as before we have, 

d 



A 



d d ke -\ (k-k ) 



X+s"k-3 +<*(qK_HUk+3 



(42) 



Exactly at q = 0, the two terms cancel so that the expectation value vanishes. For q ^ on the other hand we may 
drop the second term and write (d = 2) as before 

2 



(4f 



-e 4 



d 2 fce- A " (k - ko) " iS k+a u k _s 



1 A 2 q 2 



1 J [ i _ Vpqx 

2 4A 



d 2 ke- 



dk T e 



4A 



2A 

2 



vf ( k x - ?f 
~2Ao 



J F T2 



4Ag x 



8A 2 X 2 



e 4 



(43) 
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Transforming to position space, we have 



(27T) 
1 

8^ 2 



1 



2ttA 2 
1 

2ttA 2 



1 + i 



VFid x 
4A 



2A n A 2 





v 4A 






8A 2 A 2 







8A 2 A 2 



From this we find the phase of the order parameter has the form (recall Aq = e l9 ) 



-2 tan 



2A A 2 



(44) 



(45) 



which gives rise to a current in the — x direction, opposite to the quasiparticle current. Upon iteration of the BdG 
equations, the gap function textures in an attempt to reestablish charge conservation. 



Continuity-consistent solutions: quasiparticle wavepackets with backflow 



As demonstrated in the text, the total current has a dipolar form, 
where Vij = \ 5ij — -^r-^J is the transverse projector. 

To determine the current flow pattern in position space, we simply Fourier transform, 



J T (r) 



d?q 
(2tt) 



jT -iqr 

2 J q L 



VF 
(2tt) 2 



vf 



2tt rOO i 2 , 

dB / dq-e — ^e- 4qr 

Jo q 

2 



[Ml - yqxq v ] 



d 



d 2 



(2^) 2 



dy J y dxdy 



, 1 *V 
aq —e 4 

o q 



2tt 



dOe- lqrcosl 



d_ 

dy 



' Jo(qr) 
dq e 4 



-z x V 
-z x V 



r 



Up 

2tt 
2tt 

[x9 y - yd x ] (px(r) 



2 2 



dq J±(qr)e 



where we defined yx(r) = f^p- 
center of the wavepacket. 



(46) 



(47) 



The flow pattern is solenoidal, and falls off as 1/ r 2 far away from the 



Long-Range Current Flow Pattern from the Hall Conductance 



We have observed (see main text) that in the QH case the CS electrodynamics gives us a density of dipolar form, 



8p(r) 



1 V p -2r/ 



B . Using the multipolc expansion for the Coulomb interaction v(r) — e /r, we have for the electric 



potential due to the quasiparticle wavepacket 

V(r,6) = J rfr'<5p(r')^-^=E^r / (Mr'))^(cos(0 - 6')) 



(48) 
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Clearly by symmetry the monopole contribution (the I — term) vanishes. The dipole (Z = 1) term is nonzero, and 
thus the leading term in V is (using coordinates where angles arc measured with respect to the x-axis) 

poo p2lT 

V(r, 9) « -J-, / r' 2 dr'e- 2r ' 1b \ d0' sin(0')(cos(0 - ff)) 



(49) 

This yields an electric field of the form (valid asymptotically far from the wave packet center) 

E = —VV(r, 9) oc ^2 sin Or — cos 96^j (50) 

corresponding to a dipole moment of 7re£s/4 (note that in standard references on electrostatics, this result is often 
quoted with angles referred to the dipole axis, which is perpendicular to the choice here.) From this, using the fact 

2 

that we are in a quantized Hall state with a xy = |r we find a current flow pattern at long distances that takes the 
form 



e 


y>OC 


£ 2 B r 2 
e£ B 


Jo 


4r 2 J 


1 1 






sin 6 





„2 „3« ( 2 sin 90 -)- cos 9r 

Thus, there is an asymptotic power law tail in the current density due to the long-range nature of the Coulomb 
interaction (see also (20|.) The important points to note are the power-law form of the current and its dipolar nature, 
rather than the numerical prefactor, which is not precise since we use only the most naive estimate of the bogolon 
charge density. Further refinements will produce smaller, power law corrections to the charge density profile as well. 

I 



